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ABSTRACT 

We present new CO (1 — > 0) and (2 — > 1) observations of the N83/N84 molec- 
ular cloud complex in the south-east Wing of the Small Magellanic Cloud (SMC). 
While the (2 — > 1)/(1 — > 0) integrated line brightness ratio (in temperature units) 
is uniformly 0.9 throughout most of the complex, we find two distinct regions 
with unusually high ratios (2 — > 1)/(1 — > 0)> 2. These regions are associated 
with the N84D nebula and with the inside of the 50 pc expanding molecular 
shell N83. This shell is spatially coincident with the NGC 456 stellar association 
and the HFPK2000-448 radio continuum/X-ray source tentatively classified as 
a supernova remnant. We explore possible causes for the high ratios observed 
and conclude that the CO emission probably arises from an ensemble of small 
(R ~ 0.1 pc), warm (T g ~ 40 K) clumps. Analysis of the CO shell parameters 
suggests that it is wind-driven and has an age of slightly more than 2 million 
years. We have also used this dataset to determine the CO-to-H 2 conversion 
factor in the SMC, an especially interesting measurement because of the low 
metallicity of this source (~ 1/9 solar). Surprisingly, after comparing the CO 
luminosities of clouds in N83/N84 with their virial masses, we find a CO-to-H-2 
conversion factor Xco only 1.9 times larger than what we obtain when apply- 
ing the same algorithm to solar metallicity clouds in the Milky Way and M 33. 
This result fits into the emerging pattern that CO observations with high linear 
resolution suggest nearly Galactic values of Xco i n a wide range of environments. 



1 bolatto@astro. berkeley.edu 

2 Radio Astronomy Laboratory, Department of Astronomy, University of California at Berkeley, 601 Camp- 
bell Hall, Berkeley, CA 94720-3411 

3 Sterrewacht Leiden, P.O. Box 9513, NL 2300-RA Leiden, The Netherlands 

institute for Astrophysical Research, Department of Astronomy, Boston University, 725 Commonwealth 
Ave., Boston MA 02215 



- 2 - 



1. Introduction 

The Magellanic Clouds — with their unobscured lines of sight, their small internal ex- 
tinction, and their profusion of star-forming sites — constitute excellent laboratories for the 
study of the interaction between massive stars and their environment. They are also the near- 
est metal-poor systems with active star-formation ( Z LMC ~ Z Q /4 and Z SMC ~ Z Q /9 for 
the Large and the Small Cloud respectively; Dufour 1984; Kurt & Dufour 1998). As such, 
they provide invaluable insight into the physics and chemistry of the interstellar medium 
(ISM) in primitive galaxies. 

In this work we will discuss some of the molecular gas properties and kinematics in the 
N83/N84 region of the Small Magellanic Cloud (Henize 1956, Figure 2). This region is 
especially interesting because it is one of the few isolated, yet relatively active star-forming 
regions in the otherwise inconspicuous SMC Wing. It comprises several CO clouds (Mizuno et 
al. 2001; Israel et al. 2003), H II regions (Henize 1956; Davies et al. 1976), various early-type 
stars (Testor & Lortet 1987), and an expanding shell coincident with both the NGC 456 
stellar association and a possible supernova remnant (source 448 in Haberl et al. 2000, 
henceforth HFPK2000-448). This expanding shell appears to be interacting with the nearby 
molecular cloud, sweeping and compressing an important mass of molecular hydrogen (at 
least 10 4 M ; c.f., §3.3.2), and perhaps setting up a second generation of star formation. If 
N 83 indeed contains a supernova remnant, it represent one of a handful of opportunities to 
study the interaction between a molecular cloud and a SNR in a metal poor environment. 

Recently, Bolatto et al. (2000) have reported the detection of unusually high 12 CO 
(2 — > 1)/ (1 — > 0)~ 3 ratios in the molecular envelope of the N159/N160 complex in the 
Large Magellanic Cloud (LMC). Such high ratios are rarely observed in the Galaxy and are 
interesting because they convey information on the structure of the ISM in these systems 
(of., §3.1). Bolatto et al. suggested that these high ratios are a direct consequence of the 
low metallicity of the medium. If this is the case, we expect them to be more pronounced 
in regions with strong UV radiation fields and even lower metallicity. We undertook the 
present study in order to search for anomalous 12 CO (2 — > 1)/(1 — > 0) ratios in a region of 
very low metallicity, and thus to verify how widespread such ratios are. A second motivation 
was to attempt to further distinguish between the mechanisms that can cause such ratios. 

The present study also allows us to revisit the problem of estimating molecular masses 
using CO observations in metal-poor systems. Low heavy-element abundances cause lower 
dust-to-gas ratios in a low metallicity ISM (D/G ~ Z 1 ^ 2 ; Lisenfeld & Ferrara 1998). Be- 
cause dust particles are efficient absorbers of UV light, they shield molecules from dissociative 
radiation. Thus, the lower dust-to-gas ratios in the SMC (D/G S mc ~ 1/30 D/G M w\ Sta- 
nimirovic et al. 2000, and references therein) lead us to expect higher rates of molecular 
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photodissociation. This is particularly true for molecules such as CO, where lower metallic- 
ities also imply lower column densities, hence significantly smaller self-shielding. Indeed, it 
has long been known that emission from molecules is much fainter in dwarf galaxies, which 
tend to be low metallicity systems (e.g., Taylor, Kobulnicky, & Skillman 1998). 

Because of these effects, exactly how CO emission traces the molecular ISM in metal- 
poor systems remains uncertain, and it has stayed a goal of extragalactic CO observations 
to characterize this relationship. Since the early days of CO observations in the Galaxy 
astronomers have derived CO-H 2 calibrations and used them to compute H 2 column densities 
and masses (e.g., Sanders, Solomon, & Scoville 1984; Bloemen et al. 1986; Dame, Hartmann, 
& Thaddeus 2001, and references therein). The factor to convert 12 CO integrated intensities 
into H 2 column densities is referred to as X co = N(H 2 )/I( 12 CO), and has been determined 
using the virial theorem, X-ray shadowing experiments, and an array of other techniques. 
As it is notoriously difficult to incorporate all of the relevant physics (such as cloud structure 
and optical depth) into model calculations, astronomers have relied on empirical calibrations 
of Xco rather than on a-priori determinations. 

Since its introduction, it has been clear that X co must be a function of the local con- 
ditions of the ISM or at least the average properties of the cloud ensemble; some of the 
most relevant quantities are thought to be radiation field, dust-to-gas ratio, and metallicity. 
There are several calibrations of Xco with metallicity available in the literature (Dettmar 
& Heithausen 1989; Wilson 1995; Arimoto, Sofue, & Tsujimoto 1996; Israel 1997, 2000; Bo- 
latto, Jackson, & Ingalls 1999; Barone et al. 2000; Boselli, Lequeux, & Gavazzi 2002). These 
all agree that decreases in ISM metallicity and dust-to-gas ratio should result in noticeable 
increases in Xco, meaning a lower CO integrated intensity per unit H 2 column density. They 
differ, however, in their estimates just how large these increases in X C o will be. In partic- 
ular, observational studies find that X C o appears to depend strongly on the attained linear 
resolution. Studies based on observations with high linear resolution, such as that provided 
by interferometers, tend to find smaller increases in Xco i n metal-poor galaxies than those 
accomplished with larger beams (e.g. Wilson 1995; Walter et al. 2001, 2002; Rosolowsky et 
al. 2003). Based on comparisons with far infrared dust emission, Israel (2000) has ascribed 
these effects to intrinsic observational bias. 



2. Observations 

Carbon monoxide emission from the N83/N84 region was first detected and subsequently 
mapped by Israel et al. (1993, 2003) as part of the ESO-SEST Key Programme "CO in the 
Magellanic Clouds." Because of considerable improvements in millimeter receiver sensitivity 
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in the years since, it has become possible to use SEST 5 to make large maps of relatively faint 
sources in a reasonable period of time. It has also become possible to simultaneously observe 
the (1 — > 0) and (2 — > 1) transitions of CO, allowing us to obtain reliable line ratios. We 
mapped a region of ~ 10' x 9' encompassing N83, N84, and HFPK2000-448, simultaneously 
at 115.271204 and 230.542408 GHz (Figure 1). We detected over a dozen clouds in this area. 
The properties of the peak (2 — > 1) spectra from these clouds are compiled in (Table 1). 

We acquired the maps in November 2000, by using SEST 115/230 GHz dual SIS receiver 
and feeding the two frequencies into the split high-resolution acousto-optical spectrometer 
(AOS) backend. The resulting bandpasses (~ 50 km s _1 for 230 GHz and ~ 100 km s" 1 
for 115 GHz) easily accomodated the narrow lines and small velocity variations inside the 
molecular complex. The spectra were obtained using double-beam-switching mode, with 
a chop of 11.5' in azimuth. The few spectra contaminated by emission in the off posi- 
tion were reobserved at a different hour angle. According to its documentation, SEST has 
HPBW(230)^ 23" and HPBW(115)^ 45", and its main beam efficiencies are r/ mfe (230) « 0.50 
and r/ m b(115) ~ 0.70. We initially mapped a 24" grid (i.e., approximately full-beam spacing 
at 230 GHz) and refined this to a 12" grid in regions of bright emission and those regions 
in which we wished to increase the signal-to-noise ratio. We used a nearby SiO maser to 
determine pointing and focus corrections before each observing session. The SEST pointing 
accuracy is 3" RMS. Assuming a 63 kpc distance to the SMC, the scale of these observations 
is ~ 18 pc per arcminute. 

To compare the (2 — > 1) and (1 — > 0) data, we convolved the spectra of both transitions 
to a common angular resolution, thus removing the effect of the different beam size of the 
observations. We found unusually high 12 CO(2 — > 1)/(1 — > 0) line ratios in two regions of 
the map (Figure 2). The first instance occurs just south of N84, close to the IRAS 01133- 
7336/N84D H II region (Figure 4a). The second region is found in the gas inside the ex- 
panding shell N 83, which is associated with NGC 456 and HFPK2000-448 (Figure 4b). To 
determine the optical depth T21 of the (2 — > 1) emission we obtained 13 CO (2 — > 1) spectra 
toward several positions (Figure 5). Table 2 summarizes our observations, as well as our 
results for r 21 obtained assuming N( 12 CO)/N( 13 CO)~ 40 - 90 (c.f., §3.1.3). 



5 The Swedish-ESO Submillimetre Telescope (SEST) is operated jointly by the European Southern Ob- 
servatory (ESO) and the Swedish Science Research Council (NFR). 
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3. Discussion 

3.1. The Origin of the Unusually High Line Ratios 

Optically thick CO emission arising from thermalized gas at T > 10 K yields a brightness- 
temperature ratio (2 — > 1)/(1 — > 0)~ 1. Deviations from these conditions — for example, 
lower temperatures or subthermal excitation in gas with volume densities below the critical 
density of the (2 — >• 1) transition, n cr u ~ 10 4 cm -3 — generally produce ratios less than unity. 
However, somewhat higher ratios, (~ 1.2) are frequently found in LMC and SMC molecular 
clouds (Israel et al. 2003, and references therein). These are probably caused by a mixture 
of high and low optical depths in molecular clouds and cloud envelopes (i.e., a departure 
from the assumption of a homogeneous source characterized by a single temperature and a 
single column density). 

Unusually high 12 CO (2 — > 1)/(1 — > 0) ratios were observed by Bolatto et al. (2000) in 
the LMC, who discussed four possible origins for them: 1) self-absorbed 12 CO (1 — > 0)emission, 
2) optically thin emission from isothermal gas, 3) ensembles of small, optically thick isother- 
mal clumps, and 4) optically thick emission with temperature gradients. With the available 
data they were not able to unequivocally determine the cause of the high ratio, although 
self-absorbed CO emission appeared unlikely. The same causes may explain the high ra- 
tios we find in the N83/N84 region of the SMC. In this Section, we describe each of these 
mechanisms and explore how well they can account for the new data on the SMC. 

3.1.1. Self-absorbed emission 

Self-absorption is caused by an intervening envelope of colder CO in the line of sight, 
which produces a dip in the background CO spectra. This effect is more prominent for lower 
J transitions, since most of the intervening cold gas will be found near the ground state. 
The spectra obtained toward the region south of N84 clearly show that the 12 CO (1 -»• 0) 
transition is not suppressed by self-absorption. Fig. 4 demonstrates that the two transitions 
have the same line profile, thus convincingly showing they arise from the same parcels of 
gas. 

3.1.2. Optically thin emission 

It can be shown that for gas in local thermodynamic equilibrium (LTE) the ratio of 
optical depths in the (1 — > 0) and (2 — > 1) transitions is 
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where t\q and r 2 i are the opacities of the respective transitions, v 10 is the frequency of the 
lower transition, T ex is the excitation temperature, and k and h are Boltzmann's and Planck's 
constants. Since the opacity of rotational transitions is proportional to the square of their 
rotation quantum number J at high temperatures, optically thin warm gas can produce 
(2 — >■ 1)/(1 — »• 0) ratios < 4. The emission from such regions will be faint, and observations 
in the Galactic plane may be biased against finding them, which may explain why it is easier 
to find them in the Magellanic Clouds. However, we clearly detected emission from 13 CO in 
the same place south of N 84 where we measured a high 12 CO (2 — > 1)/(1 — > 0)= 2.0 ± 0.2 
ratio (Figure 5). Although the observed 12 CO/ 13 CO (2 — ► 20 ratio is somewhat high, it 
nevertheless indicates that this transition is optically thick unless the carbon isotopic abun- 
dance ratio is much lower in the SMC than in the Galaxy. The best available determinations 
of isotopic abundance ratios in the SMC suggest that N( 13 C)/N( 12 C)~ 40 — 90, although 
data are scarce (Heikkila, Johansson, & Olofsson 1999). In Table 2 we summarize our op- 
tical depth results for several clouds in the complex. The range of r 2 i values listed in the 
last column are calculated assuming the quoted limits from Heikkila, Johansson, & Olofsson 
(1999). 

For the pointing at (+0.0, —0.4) we find an opacity in the range r 2 i ~ 1.5—4.4, indicating 
that the (2 — > 1) emission there is optically thick. If the N( 12 CO)/N( 13 CO) column density 
ratio turns out to be significantly lower than 40, however, the emission could be optically thin. 
Such a scenario could conceivably result from the combined effects of chemical fractionation 
and selective photodissociation in the ISM. Heikkila, Johansson, & Olofsson (1999), for 
example, determined this ratio to be ~ 15 in the N27 region of the SMC using a mean 
escape probability radiative transfer code modeling spherical clouds. 



3.1.3. Emission from very small warm clumps 

Collections of very small, isothermal, optically thick clumplets also have the potential 
to produce high 12 CO (2 — > 1)/(1 — > 0) ratios. They do so partially through a beam filling 
effect; because of the faster growth of the (2 — > 1) opacity with column density the optically 
thick region of the clumps appears larger in the higher transitions, thus filling more of the 
beam. Figure 6 illustrates the geometry of this effect, in which the optically thick r 2 i = 1 
"photosphere" of the clump is larger and occurs closer to its surface than the t w = 1 
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"photosphere." Perhaps more importantly, small clumps provide a natural geometry that 
allows the (2 — > 1) transition to become optically thicker than the (1 — > 0) transition. The 
line ratios for a spherical geometry can be easily computed under the assumption of local 
thermodynamic equilibrium (LTE); this is probably a safe assumption if the clumps have 
volume densities in excess of n > 10 4 cm -3 (because of the critical density n cr i t oc J 3 , it 
is extremely difficult to obtain a high (2 — > 1)/(1 — > 0) ratio with subthermal excitation 
unless the equivalent radiation temperature is higher than the gas temperature, which is not 
usually the case in mm-wave radio astronomy). For each transition, we assumed a constant 
gas temperature and density and numerically integrated to find the opacity, r(r), in a clump 
as a function of r, the projected distance to the center. The results of these calculations 
as a function of the gas temperature, T g , and the central clump opacity, t° , are shown in 
Figure 7. The shaded regions represent the observational constraints: the measured 12 CO 
(2 — > 1)/(1 — > 0) ratio and the central clump opacity derived from the 12 CO/ 13 CO (2 — > 1) 
measurements as described in the previous paragraph. The observations constrain the gas 
temperature to be strictly T g > 30 K and more likely T g > 40 K. The clump radius is 
constrained to be r c ~ 1 in opacity units for the (2 — > 1) transition. This opacity translates 
into a physical column density iV( 12 CO) ~ 1.5 x 10 17 cm~ 2 at T g ~ 40 K, which can be 
converted to a length of R ~ 0.1 pc by assuming that all the carbon is locked in the CO 
molecules with a carbon abundance of log(C/H) k, —4.7 (Garnett et al. 1995; Kurt & Dufour 
1998) and that the volume density is nn 2 ~ 10 4 cm -3 (i.e., about the critical density of the 
(2 — >• 1) transition). The resulting parameters thus indicate fairly standard, but very warm, 
clumps. The heating is presumably provided by the star cluster found near the high-ratio 
region (source 166 of Bica & Schmitt 1995) or by the adjacent H II region and IR source 
(IRAS 01133-7336; Figure 2). 

The appeal of the "small warm clumps" scenario is twofold: 1) it naturally explains 
why it is easier to find high line ratios in systems of low metallicity, such as the LMC and 
SMC, and 2) as it has predictive power, it can easily be falsified. Indeed, we have shown that 
clump radii have to be small in opacity units to produce large ratios. If the physical sizes 
of clumps do not change with metallicity, low CO column densities will be more common 
in lower metallicity systems because of the lower abundances of carbon and oxygen. If the 
abundance of CO scales by a factor similar to that of the metallicity, then a clump with 
central opacity r c ( 12 CO) ~ 4 in the Milky Way would have r c ( 12 CO) ~ 0.5 in the SMC. 
The "small warm clumps" model requires that the 13 CO (2 — > 1)/(1 — > 0) ratio be ~ 3 — 4 
in this source, and it can be used to predict the intensities of the higher CO transitions. The 
model also requires that the N( 12 CO)/N( 13 CO) column density ratio be ~ 40 — 50 south of 
N84. 
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3.1.4- Optically thick regions with large temperature gradients 

High brightness ratios can, in principle, also be produced by optically thick gas if the 
t = 1 surfaces of each transition occur in regions of different temperature. To reproduce 
the observations, we need a ratio of excitation temperatures T ex (r 2 i) /T ex (T W ) > 2. How 
large a temperature gradient does this change represent? Using the LTE approximation 
we computed the column densities of CO necessary to attain r = 1 for each transition at 
different temperatures. We used this information to obtain a very rough estimate of the 
typical increment in column density AiV( 12 CO) required to go from the T21 = 1 surface at 
T, to the Tio = 1 surface at T/2. We found that a reasonable estimate is AiV( 12 CO) ~ 10 16 
cm -2 . Taking into account the abundance of carbon, this is equivalent to a change in visual 
extinction AAy ~ 0.05. Since CO becomes the main reservoir of carbon at Ay > 1, this 
means that the cloud turns optically thick in the (1 — > 0) transition at a depth only ~ 5% 
larger than that at which the r 2 i = 1 surface is found. To explain the observations, however, 
we require the excitation temperature of the gas to drop by a factor of two for this ~ 5% 
change in depth. Gas temperature profiles computed for PDRs externally heated by UV 
sources show that such gradient is extremely unlikely (e.g., Hollenbach & Tielens 1999). 
This simple argument is also supported by line ratios obtained from self-consistent PDR 
calculations. At no point in the parameter space of a homogeneous, externally heated PDR is 
the CO (2 — > 1)/(1 — >■ 0) > 1.8, and ratios greater than unity only occur at very high densities 
and radiation fields (Kaufman, Wolfire, Hollenbach, & Luhman 1999). Perhaps embedded 
IR sources with no UV flux, such as protostars, could produce the heating necessary to 
obtain the right temperature ratio without the associated photodissociation of CO. This 
appears very unlikely, however, since these sources would have to occur at a particular 
depth, coincident with the r 2 i = 1 surface. 

We have explored four possible explanations for the unusually high 12 CO (2 — > 1)/(1 — > 0) 
ratios observed south of the N84 region. We have shown that two mechanisms, self- 
absorption and optically thin emission, are improbable. Of the remaining two, small warm 
clumps and optically thick gas with large temperature gradients, we prefer the first, although 
it is not possible to entirely discard the second with the present data. High line ratios could 
also conceivably originate in the optically thin portion of clumps with density gradients (for 
example, critical Bonnor-Ebert spheres). The proper modeling of such objects, however, 
requires a full self-consistent radiative transfer and temperature profile calculation, and is 
beyond the scope of this paper: nevertheless we note that warm temperatures will still be 
necessary to produce high line ratios. In addition, the relevant ratios are also obtained if 
we allow completely inhomogeneous models. One of the simplest examples would be a pop- 
ulation of cold, dense clumps surrounded by a warm, tenuous envelope. Unfortunately, a 
model containing only two completely independent components already has a minimum of 
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eight free parameters, which greatly exceeds the number of constraining observational pa- 
rameters (three in our case). We may decrease the number of free parameters somewhat by 
making plausible assumptions, but even in that case would we need more constraints than 
we presently have. Future observations may remedy this. 

3.2. The Expanding Molecular Shell 

The second instance of unusually high 12 CO (2 -> 1)/(1 -»• 0) ratios in the N83/N84 re- 
gion occurs in a remarkable place: the inside of the expanding shell-like structure associated 
with NGC 456 and HFPK2000-448 (Figure 2). This optical shell has a molecular coun- 
terpart apparent in a careful inspection of Figure 3. Starting at velocities 155-157 km s _1 
the center of N83 lights up, while at 161 km s -1 the eastern rim emits brightly in CO. At 
165 km s _1 , some emission shows up at the center of the shell in the maps of the (2 — >■ 1) 
transition. The spectrum obtained over the central 1' of this structure even more clearly 
exhibits this kinematic signature (Figure 4b). While the emission associated with the walls 
of the shell has very similar brightness temperature in the (1 — > 0) and (2 — > 1) transitions, 
the gas inside the shell appears brighter in the (2 — >• 1) transition by a factor of ~ 2. The 
shell is ~ 50 pc in diameter, and the molecular material associated with it is expanding at 
a velocity ~ ±6.5 km s -1 . The northern and eastern edges of the shell are sharp and well 
defined in the optical pictures, while its southwestern edge appears ragged and soft. The 
cause of this difference in appearance is strikingly clear in the CO channel maps: along the 
northern and eastern edges, the expansion of the shell is contained by a molecular cloud at 
Vlsr ~ 163 km s -1 , while the south and west sides appear able to expand more freely. It 
is not immediately clear whether this bright CO cloud (corresponding to peaks MP 5 and 
MP 7 in Table 1) consists of swept-up matter, or represents an ambient cloud interacting 
with the expanding shell. 

A question arises: is the expanding shell due to a supernova explosion, or is it a wind- 
driven structure? The classification by Haberl et al. (2000) of HFPK2000-448 as a SNR 
was tentative, and based on the detection of coincident radio continuum emission at 13 cm 
and X-Ray emission by ROSAT. Extended X-ray emission and/or extended nonthermal 
radio continuum are telltale evidence for SNRs. However, in the case of HFPK2000-448 
the data are not conclusive. The X-ray emission could be extended but the quality of this 
determination (as measured by its likelihood parameter) is low, and there is no detection of 
the radio continuum at other frequencies to confirm that the spectral index is nonthermal. 

In any case, the rather modest expansion measured in CO is hard to reconcile with 
the much higher expansion speeds expected from a supernova remnant. It is quite in line, 
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however, with the much lower speeds that should characterize wind-driven shells. From the 
interstellar bubble model by Weaver et al. (1977), we find that the age, t, of the bubble 
is related to its radius, R s , in parsecs and its expansion velocity, V s , in km s _1 by t = 
6 x 10 5 i? s /Vs. With a radius of about 25 pc and an expansion velocity of 6.5 km s _1 , the age 
of the bubble should be about 2.3 million years. From the model we also find the relation 
between mechanical wind-luminosity, L w , in erg s" 1 and original ambient total hydrogen 
density, n H) in cm -3 to be L w /n H = 7.65 x 10 46 i? 5 /t 3 . By conservatively assuming L w / L hol = 
10~ 3 for the ratio between mechanical wind and stellar bolometric luminosities, we may 
reduce this to L bo i/n H = 6 x 10 37 ergs _1 cm 3 = 1.5 x 10 4 L o cm 3 . If the total stellar luminosity 
of the embedded association NGC 456 is of the order of 10 7 L , the original ambient density 
should have been of the order nn ~ 650 cm' 3 . These luminosities and original densities seem 
quite reasonable for a luminous association dispersing its molecular birthcloud. It is also not 
unreasonable to suppose that a massive member of the association turned supernova within 
a few million years. If such an event indeed occurred, the resulting remnant should rapidly 
have filled the wind-blown bubble that we have observed. Consequently, the chemistry 
and excitation around N 83 could be dominated by SNR-caused shocks, even though the 
overall kinematics of the shell would still be those of a wind-driven structure. We found 
no evidence, however, of uncommon mm-line ratios in the N83B/C molecular cloud. Its 
12 CO (2 — > 1)/(1 — > 0) ratio is ~ 0.9, identical to that of the rest of the complex, and our 
observations of HCO+ (1 -> 0) toward offsets (-3.6,4.6) yield a ratio 12 CO/HCO+ « 18 ±2, 
similar to that reported by Heikkila, Johansson, & Olofsson (1999) toward the SMC source 
N27. 

It would be of considerable interest to confirm the SNR identification of N 83 /HFPK2000- 
448. If it is confirmed, N 83 in the SMC would be very similar to the well-established SNR 
N 49 in the LMC, which also appears to be interacting directly with a molecular cloud (Ba- 
nas et al. 1997). N 83 would represent another excellent opportunity to study molecular 
cloud-SNR interactions in a metal-poor environment, complementing studies of similar in- 
teractions within the Galaxy. Indeed, notwithstanding the fact that most Type II supernovae 
must occur within their parent clouds and despite their potential importance as triggers of 
a secondary wave of star-formation activity (Opik 1953), there are only a handful of such 
molecular cloud-SNR associations identified in the literature (e.g., Cornett, Chin, & Knapp 
1977; Wootten 1977; Wilner, Reynolds, & Moffett 1998; Kim et al. 1998; Dubner et al. 1999, 
and references therein). Most of these interactions lie in the Galactic plane, where obscura- 
tion and confusion pose observational difficulties that are largely absent in the Magellanic 
cluds. 
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3.3. The X co Factor in N83/N84 

3.3.1. Measuring Cloud Properties 

To measure the individual cloud properties we proceeded as follows: from the individual 
map spectra we constructed uniformly sampled main beam intensity data cubes in the CO 
(2 — > 1) and (1 — > 0) transitions using the COMB built-in 2D gaussian interpolation. The 
cubes have angular resolutions (HPBW) of 35" and 55", respectively, corresponding to spatial 
resolutions of 11 pc and 17 pc. Their velocity resolution is 0.25 km s _1 . After convolving 
the (2 — > 1) cube to the resolution of the (1 — > 0) cube, we carried out a pixel-by-pixel 
comparison of the two cubes and found that the best linear relation between both datasets 
was Ico(l ~~ ► 0) = 0.9 Ico(2 — > 1), very similar to other results obtained for Magellanic Cloud 
objects (c.f. Israel et al. 2003, and references therein). This scaling will be used later to 
compute molecular masses using CO (2 — >• 1) observations. 

We proceeded to identify clouds in the channel maps, and then measured three properties 
for each cloud: its size, its velocity dispersion, and its integrated CO intensity. From the 
first two quantities we calculated a virial mass for each cloud, which we compared to the 
mass derived from the integrated intensity using the Galactic CO-to-H 2 conversion factor. 
To obtain these numbers we specified a threshold intensity, then built a mask that started 
from a seed position and velocity, and included neighboring regions of the data cube with 
two adjacent velocity channels above this threshold. We defined a "cloud" as an unmasked 
region contiguous in a — 5 — v space. This is a practical definition that is similar or identical 
to those employed in a number of other studies (e.g. Solomon et al. 1987; Heyer et al. 2001; 
Rosolowsky et al. 2003). The threshold intensity was varied over a wide range, and cloud 
properties were computed for each of its values. Obviously, very high values of the threshold 
will select only small regions within a cloud, while low values will tend to merge clouds 
together, and even lower values will include positive noise as emission. We determined the 
reasonable range of values to use for each cloud by observing when the cloud properties 
experienced a "jump" caused by the mask growing into a neighboring cloud. 

After identifying an individual cloud, we calculate its size, velocity dispersion, and 
CO flux as follows. To obtain the size, we compute the second moment of the integrated 
brightness distribution of the cloud in the a — 5 plane (i.e., the RMS size). Because the 
moment underestimates the true size for gaussian clouds, we include a correction factor that 
depends on the ratio of the threshold intensity used to the peak intensity of the cloud. The 
inclusion of this factor is important for two reasons: 1) in cases of confusion and low signal- 
to-noise, it allows us to estimate the properties of an entire cloud from the portion that is 
easily identifiable, and 2) it avoids oversubtraction in our beam deconvolution. It assumes, 
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however, that the cloud is gaussian in position-velocity space. Fortunately, this assumption 
can be veryfied by observing how the cloud properties change as a function of threshold 
intensity (they should be constant), and it appears to hold for the clouds here studied. From 
the corrected second-moment values, we "deconvolve" the beam by subtracting its size from 
the measured cloud size in quadrature. Finally, we apply the factor of ^7= (Solomon et al. 
1987) to convert from the RMS size to the radius of a theoretical spherical cloud (this factor 
can be readily motivated by considering a constant density spherical cloud and comparing 
the RMS size to the spherical radius). Along with an assumed SMC distance of 63 kpc (see 
Walker 1998; Cioni et al. 2000), this treatment results in a physical cloud size given by the 
formula: 



R sph * OAly/al + al-2al am , (2) 

where a a and as are the threshold-corrected RMS sizes of the cloud in the a and S directions, 
&beam is the RMS size of the beam (all measured in arcseconds), and R sp h is the spherical 
radius of the cloud in parsecs. 

Other works, in particular those focused on Milky Way clouds, have treated clouds as 
ellipsoids and used the corresponding form of the virial theorem due to Bertoldi & McKee 
(1992). This approach is problematic in extragalactic astronomy, because the observations 
do not have many resolution elements across the cloud and their signal-to-noise is low. 
Under these conditions deconvolution becomes very difficult and the precise shape of the 
source is hard, if not impossible, to recover. Therefore, we content ourselves with measuring 
a single size parameter and considering the simple model of a spherical cloud. 

We calculate both the threshold-corrected RMS velocity dispersion, a V) and the equiv- 
alent width, EW, from the integrated spectrum of each cloud. For the equivalent width, we 
employ the definition of Heyer et al. (2001), 



EW = ^ {Vi) (3) 
Max(VK)) { ' 

where ipivi) is the spectrum of the cloud in question. For a gaussian distribution, 
EW = 1.06 FWHM. Because it does not diverge when signal-free noise is added to the 
edge of a spectrum, this linewidth measure is more robust than the second moment as a 
measure of the RMS velocity dispersion. Moreover, by comparing these two measures we get 
a good idea of how well a spectrum is approximated by a gaussian (for our data, the ratio 
EW/cr v correlates very well with the y 2 value of a gaussian fit to a spectrum). 



We compute the CO luminosity by summing the intensities of the entire cloud and 
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again applying a correction factor based on the peak-to-threshold ratio and an assumed 
gaussian profile. In the case of the CO (2 — > 1) data cube, we applied the aforementioned 
0.9 conversion factor to calculate the corresponding CO (1 — > 0) luminosity for the cloud: 
essentially, we used the (2 — > 1) measurements as predictors of the (1 — > 0) intensities, taking 
advantage of their higher S/N and spatial resolution. To compute molecular cloud masses 
(including the contribution by He) we use the formula 

M mo i = a C o S C od 2 , (4) 

where Sco = / Ico dfl is the spatially integrated 12 CO flux of the molecular cloud measured 
in Jy km s _1 , a C o — 1-07 x 10 4 M Mpc~ 2 (Jy km s -1 ) -1 (corresponding to a canonical 
conversion factor X co = 2 x 10 20 cm" 2 (K km s -1 ) -1 ), and d is the distance to the SMC in 
Mpc. 

We compute the virial masses using the assumption that each cloud is spherical and 
virialized with a density profile of the form p oc r -1 . Thus the virial mass is given by the 
formula (Solomon et al. 1987) 

M vir = 1040 a 2 v R sph , (5) 

where a v is the RMS cloud velocity dispersion in km s -1 , R sp h is the spherical cloud radius 
in pc, and M vir is the cloud's virial mass in M Q . 

3.3.2. Results 

We have applied the procedure described above to several subsets of the data, as well 
as the whole datacube. We focused our attention on three clouds contained in the Nyquist- 
sampled regions of our original dataset, roughly corresponding to MP7, MP10, and MP13 in 
Table 1 (i.e., the clouds associated with N83B/C, N84C, and N84B/D respectively; Fig. 2). 
We also computed masses for what we call the "northern region" (the complex encompassing 
N83B/C and N84C but not including MP6, which appears to be an unrelated cloud peaking 
a much lower velocities), and for the complex as a whole. Our results are summarized in 
Table 3. Because the dominant sources of uncertainty in our determination of cloud masses 
are the cloud boundary definitions, we repeated the procedure described above for a wide 
variety of thresholds. Thus, we probed the properties of the complex on a variety of size 
scales, and on each scale we estimated a reasonable range of values for each of the physical 
parameters. In Table 3 we give the properties of structures on scales ranging from what 
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seem to be compact, self-gravitating object to complexes composed of a few large, bound, 
virialized clouds connected by low surface brightness filaments. 

How do the results obtained from the (2 — > 1) data compare with those from the more 
commonly observed (1 — > 0) transition? CO (1 — > 0) usually provides the basis for estimating 
X C o in the SMC and other galaxies, even though the lower spatial resolution at 115 GHz 
makes size measurements more difficult. We performed identical analyses on both datasets. 
For the most part, the 115 GHz properties were similar to those derived from the 230 GHz 
data, though the CO flux and linewidth of MP10 are both notably larger at 115 GHz. The 
complex, as a whole, has somewhat higher linewidth (and consequently virial mass) when 
measured in the (1 — > 0) transition. This is artificial and caused by the higher noise of the 
(1 — >• 0) datacube: to encompass the same area used in the (2 — > 1) calculations we need to 
use too low a threshold, causing positive noise to be included as emission. 

Because of the higher resolution and better signal-to-noise, we prefer the results from 
the (2 — -> 1) data. Assuming all clouds share the same conversion factor, these results yield 
X co « 7.0 ± 3.4 x 10 20 cm" 2 (K km s" 1 )- 1 (i.e., 3.5 ± 1.7 times the Galactic value) on the 
size scale of individual molecular clouds (~ 10 pc). Using this conversion factor, the total 
molecular mass traced by CO in the complex is M(H 2 ) ~ 2.5 x 10 5 M . 

We believe our algorithm to be both robust and simple. It is designed to be easily 
applicable to a heterogeneous set of observations. However, because several details (particu- 
larly the masking correction) differ from those previous works, our results (particularly our 
size measurements) may not be readily comparable to those in the literature. To account 
for this, we applied our algorithm to both the Rosette Molecular Cloud (Williams, Blitz, & 
Stark 1995) and a subset of well resolved solar metallicity clouds from the BIMA survey of 
M 33 (Rosolowsky et al. 2003). Both of these tests yielded values for X co of 3.4 - 4.0 x 10 20 
cm~ 2 (K km s -1 ) -1 . In light of this, our result can be restated as follows: our measurments 
of the Xco factor in the SMC yield a value for Xco 1.9 ±0.9 times the value of Xco obtained 
for identical measurements on the Rosette Molecular Cloud and several M 33 clouds. 

How do these molecular masses and estimates of Xco compare to those obtained by 
other studies of the SMC? The results presented here are consistent with those obtained 
independently by Israel et al. (2003), taking into consideration that they mapped only part 
of the complex in the (1 — >• 0) CO transition, with generally lower signal-to-noise ratios. 
Our (2 — > 1) result is lower by about a factor of ~ 2 than the average value X C o~ 12 x 10 20 
cm" 2 (K km s" 1 )" 1 obtained by Rubio et al. (1993) from SEST CO (1 -»• 0) measurements 
of clouds elsewhere, in the southwest Bar of the SMC. However, when we adjust the sizes of 
Rubio et al. (1993) to match our own definition and include beam deconvolution, we find a 
~ 40% reduction in the average size of a cloud. The corresponding drop in cloud virial masses 
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suggests a value of X co ~ 6 — 8xl0 20 cm -2 (K km s -1 ) -1 consistent with our estimates. We 
also take the opposite approach, bringing our results into line with their measurements 
by considering the properties of the N83 complex at the 0.4 K km s -1 integrated intensity 
contour. If we reject the low velocity emission, as Rubio et al. (1993) did for several clouds, 
and employ their definitions and formulae then we obtain a conversion factor of Xco~ 
20 x 10 20 cm -2 (K km s -1 ) -1 , in excellent agreement with their average result for complexes 
(23 x 10 20 cur 2 (K km s -1 ) -1 ). This value is also in agreement with that obtained by Mizuno 
et al. (2001) using the NANTEN telescope, with a beam size HPBW^ 2'.6. 



3.3.3. Caveats and Limitations 

We note that most common limitations inherent to the data such as low signal-to-noise, 
blended line profiles, and low spatial resolution will drive our estimate of Xco toward higher 
rather than lower numbers. All the usual caveats regarding the applicability of the virial 
theorem apply, with the added concern that we are considering a particularly active region. 
The effects of an expanding shell and active star formation on nearby clouds, although not 
immediately apparent from the data, may be important. 

We should caution that virial estimates of X co on the large scales are suspect. First, 
they require assuming that the entire complex is virialized and that large scale motions 
(e.g., galactic rotation, expanding shells) are negligible. Second, the mass on the spatial 
scale of the measurement has to be dominated by molecular hydrogen, meaning that the 
contributions from stars and H I are negligible. The ATCA+Parkes HI observations of this 
region (Stanimirovic et al. 1999), for example, show that the atomic hydrogen mass in a 120 
pc diameter region centered on the N83/N84 complex is ~ 530 x 10 3 M , very similar to 
our virial mass estimate. Because this is an integral along the line-of-sight, however, it is 
impossible to assert how much of this mass is really associated with the complex. Third, 
they require assuming that the velocity field inside the region is uniformly traced by CO 
emission, otherwise the intensity-weighted second moment usually employed to obtain the 
velocity dispersion is meaningless. 

Finally, we have gauged the magnitude of the conversion factor indepently from virial 
considerations by following the method outlined by Israel (1997). We use his SMC calibration 
N H /a(FIR) = 1.65 x 10 28 cm- 2 (Wm- 2 sr- 1 ) 1 to derive the molecular mass from the far 
infrared emission near the N83/N84 complex. At 15' resolution, the temperature-corrected 
far-infrared surface brightness is a (FIR) = 7 x 10" 7 Wm" 2 sr- 1 (Schwering 1988; Schwering 
& Israel 1990) so that N R = 1.2 x 10 22 cm~ 2 . Both the old H I map by McGee & Newton 
(1981) and the newer ATCA+Parkes data by Stanimirovic et al. (1999) yield a mean neutral 
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hydrogen column density of A/hi = 5.5 x 10 21 cm~ 2 . From this we deduce the molecular 
hydrogen column density A/h2 = 3.2 x 10 21 cm -2 , corresponding to a molecular gas fraction 
2A/h2/(A/h + A^He) — 0.4, which is not unusual for a star-forming region. The CO luminosity 
integrated over the N 83/N 84 complex is given by Mizuno et al. (2001) as Lqo — 1-3 x 
10 4 K km s _1 pc 2 . As this luminosity is beam-independent for unresolved sources, we can 
easily extrapolate to the 15' resolution of the molecular hydrogen column density just derived. 
At that resolution, the CO integral corresponding to the observed CO luminosity would be 
-^co = 0.22K km s -1 . This is probably an underestimate, because such a large beam would 
encompass additional H II regions most likely associated with their own CO clouds. Our 
best estimate - to be verified by future CO mapping of the area - is therefore I C o — 0.33 ± 
0.11K km s -1 , yielding an independent estimate Xco — 100 ± 50 x 10 20 cm" 2 (K km s" 1 ) -1 
very similar to the other X co values of the SMC found in the same way (Israel 1997). 

These results fit well into an apparently consistent pattern. Observations with high 
linear resolutions, i.e. single-dish observations of nearby extragalactic objects such as those 
presented here or interferometric observations of more distant objects, consistently suggest 
virial cloud masses implying nearly Galactic Xco conversion factors regardless of environ- 
ment. By contrast, observations not resolving individual clouds but covering whole complexes 
have suggested much higher conversion factors in metal-poor systems. This pattern may be 
explained in several ways: 1) it may be caused by a large "mm-dark" molecular mass in 
metal-poor galaxies, comprised by H 2 clouds without CO emission. In this scenario the 
observed CO emission arises only from the densest, most shielded clouds in a complex. If 
this is true, CO observations suffer from a severe observational bias. 2) It is possible that the 
complexes observed are not virialized, or their masses dominated by unaccounted-for stellar 
or atomic gas components, and that there is not a great deal of hidden molecular mass. 
In this case, high-resolution observation such as those presented here, suggest that individ- 
ual clouds are virialized with an approximately Galactic X<x>, indicating that metal-poor 
galaxies have molecule-poor ISMs. 

If the complexes are virialized and their masses are dominated by the molecular gas, 
then their global dynamics would be better probes of the overall molecular mass than the 
dynamics of individual clouds. On the other hand, if complexes are not virialized or their 
masses not dominated by molecular gas, then there may be no way for CO observations 
to get a handle on the total molecular masses of metal-poor galaxies. In this case other 
tracers, such as far-infrared emission from dust, large surveys of UV H 2 line absorption, or 
FIR [C II] measurements, may be the only way of obtaining global molecular mass estimates. 
By themselves, CO observations are insufficient to distinguish between these two scenarios. 
Observations of other tracers, such as the far-infrared, will be necessary to establish whether 
CO is a good tracer in metal-poor galaxies. 
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4. Summary and Conclusions 

We have presented new CO (1 — > 0) and (2 — >• 1) observations of the N83/N84 nebulae in 
the SMC that show an active and complex region, home to several interesting astrophysical 
phenomena. We have identified 14 molecular clouds and an expanding molecular shell with 
a center coincident with the NGC 456 stellar association and the HFPK2000-448 radio 
continuum/X-ray source classified by Haberl et al. (2000) as a SNR. The shell is most 
likely wind-driven, and appears to have an age of 2.3 million years. We have also identified 
two extended regions of unusually high (2 — > 1)/(1 — > 0) ratio: the first found south of 
the complex near the N84D emission line region, the second located inside the expanding 
molecular shell. A 13 CO spectrum toward this first region shows that the (2 — > 1) transition 
has opacity r « 1.5 - 4.4 assuming N( 12 CO)/N( 13 CO)~ 40 - 90. 

Barring the unlikely cases that the isotopic N( 12 C)/N( 13 C) carbon ratio is consider- 
ably lower than in the Galaxy or that fractionation effects have decreased the isotopomer 
N( 12 CO)/N( 13 CO) ratio to much below the isotopic ratio, this high line ratio is not caused 
by low optical depths. Inspection of the line profiles in Figure 4 also shows that the high 
ratio is not due to self-absorbed (1 — » 0) emission. We have considered optically thick 
emission with temperature gradients, but some simple calculations show that the gradients 
required to reproduce the observed ratio are unrealistically large compared to those expected 
in photodissociation regions. Finally, we modeled the emission produced by an ensemble of 
spherical, isothermal clumps, and found that we can reproduce the observations with fairly 
warm (T g ~ 40 K) but otherwise normal molecular clumps (n ~ 10 4 cm -3 , R ~ 0.1 pc). 
The high temperature required may explain why these high ratios are not more commonly 
observed. 

We have also used the CO data to revisit the problem of the CO-H2 conversion factor 
in the SMC. Using a simple and robust algorithm we found the Xco factor in three clouds 
of the complex to be X co « 7.0 ± 3.4 x 10 20 cm~ 2 (K km s" 1 )" 1 (X co « 3.5 ± 1.7 times the 
Galactic value), estimated from the CO (2 — > 1) data using virial calculations and supported 
by the results obtained from the (1 — » 0) data. This value is about 1.9 times that which we 
obtain when applying our algorithm to the Rosette Molecular Cloud and GMCs in M 33. 
The result is an Xco that is considerably smaller than would be expected a priori for the very 
metal-poor SMC, which has a metallicity Z ~ 1/9Z & . Higher values have been obtained 
from a CO/virial analysis of the complex as a whole, and higher yet from a non- virial, far- 
infrared analysis of the complex. These results are consistent with the pattern emerging 
from extragalactic observations over the last decade. 
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Table I. CO (2 -> 1) Molecular Peaks in the N83/N84 Region 
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Table 2. Optical Depth Data 
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0.16+0.02 


160.63+0.17 


4.58+0.41 


0.82+0 


06 


7.6 


- 17.0 


-2.8 


+3.8 


12 co 


(2 


-1) 


3.04+0.02 


162.45+0.02 


3.07+0.03 


9.84+0 


14 










13 co 


(2 


1) 


0.4+0.02 


162.05+0.05 


2.37+0.11 


0.96+0 


06 


4.0 


- 9.2 


-3.6 


+4.6 


12 co 


(2 


-1) 


3.5+0.02 


161.87+0.01 


2.57+0.01 


9.34+0 


12 










13 co 


(2 


-1) 


0.4+0.02 


161.81+0.03 


2.36+0.08 


0.98+0 


02 


4.4 


- 10.0 
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Table 3. Cloud Properties 



Quantity MP 7 MP 10 MP 13 N. Region c Complex 



CO (2 -> 1) 



«off 




-3^09 


-175 


+0:02 


-3^00 


-2!34 


5 off 




+3^98 


+5^01 


+0:39 


+4^35 


+3^62 


Velocity Centroid (km s _1 ) 




162.0 


160.3 


168.3 


161.6 


162.0 


Velocity Dispersion a v (km s" 




1.25 


1.45 


1.1 


1.94 


4.7/3.0 a 


Rsph (pc) 




22 


14 


13.5 


32 


65 


S co (Jy km s- 1 ) b 




497 


136 


115 


742 


1,700 


M mol from CO (10 3 M ) 




21 


5.8 


4.9 


32 


73 


M vir Virial Mass (10 3 M ) 




35 


29 


18 


126 


1,500 


M vir /M mol 




1.7 


5.0 


3.7 


4.0 


20 


«CO d 
Ctco(MW) 




0.9 


2.7 


1.9 


2.2 


11 


CO (1 -> 0) 


«off 




— 3fl0 


-173 


+0!24 


-2:92 


-lf91 






+4^02 


+5^03 


+0:54 


+4^43 


+3^64 


Velocity Centroid (km s _1 ) 




162.2 


160.0 


168.4 


161.5 


162.2 


Velocity Dispersion a v (km s" 




1.2 


2.0 


1.25 


2.2/1. 8 a 


4.2/3.2 a 


Rsph (pc) 




21 


15 


16 


35 


72 


Sco (Jy km s _1 ) 




440 


193 


132 


715 


1,500 


M mol from CO (10 3 M ) 




19 


8.2 


5.6 


30 


63 


M vir Virial Mass (10 3 M ) 




32 


60 


25 


176 


1,300 


M vir /M mo i 




1.7 


7.3 


4.5 


5.8 


21 


«CO d 

a co (MW) 




0.9 


3.9 


2.4 


3.1 


11 



a Dispersion from the second moment/from the gaussian equivalent width. The two differ 
because the line is not well described by a single gaussian. Unless specified, the dispersion 
referred to in the text is that obtained from the second moment. 

b The 12 CO (1 — > 0) intensity is derived from the (2 — > 1) measurements, using 
Ico(l->0) = 0.9I C o(2->l). 

c The Northern Region is defined to include the clouds associated with N84C and N83B/C, 
but excludes MP6. 

d Computed by dividing the previous row by the ratio obtained using the same algorithm 
in Galactic and M 33 sources (1.85). 
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R.A.offeet 

Fig. 1. — Integrated intensity CO (2 — 1) map, showing the sampling. The contour levels 
start at 1 and go in steps of 1 to 8 K km s _1 . This map was produced using gaussian 
interpolation, smoothing to a final beam size of 38" HPBW. The offsets are with respect to 
a 1950 = l h 13 m 23U, 5 1950 = -73 36'33" (a 200 o = l h 14 m 4K3, <5 200 o = -73°20'42"). 
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Fig. 2. — CO (2 — >■ 1) integrated intensity map overlaid on the annotated DSS optical image 
of the N83/N84 region. The different symbols indicate: IR sources in the SIMBAD database 
(stars), molecular clouds listed in Table 1 (crosses), and individual pointings with high 12 CO 
(2 — >■ 1)/(1 — > 0) ratio (circles). The labels indicate the approximate positions of the emission 
line nebulosities listed in the Henize (1956) catalog, as well as two sources associated with 
the expanding molecular shell: the NGC 456 stellar association and the HFPK2000-448 radio 
continuum peak/SNR. The white circle in the bottom right corner indicates the resolution 
of the CO map, identical to that of Figure 1. 
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Fig. 3. — Channel maps of the CO emission for the (1 — > 0) (top group of twelve panels) and 
(2 — >■ 1) (bottom) transition, convolved to the same angular resolution (55") and overlaid on 
the DSS image of the region. The contours are identical for both transitions: starting at 
0.2 and increasing by 0.2 K km s _1 . The ratio of main beam temperatures is uniformly 1 
throughout the map, except in the two regions indentified in Figure 2. 
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Fig. 4. — (a) CO (1 — >■ 0) (thick line) and (2 — > 1) (thin line) toward the region south of N84 
(N84D), corrected by their respective main beam efficiencies and with the higher transition 
convolved to the angular resolution of the lower one. Remarkably similar velocity structure, 
and a line ratio ~ 2.0. (b) Same toward the CO expanding shell, with both transitions 
convolved to 1' resolution. The (2 — >■ 1) transition is enhanced in the center component, 
corresponding to gas filling the inside of the shell. 
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Fig. 5. — 12 C0 and 13 C0 (2 — > 1) spectra obtained near several peaks in the integrated in- 
tensity map. The pointing near the N84D region with offsets (+0.0, —0.4), where we observe 
12 CO (2 — > 1)/(1 — > 0) ratio ~ 2.0, was also observed. This position was chosen because of 
its relatively strong 12 CO emission. Unfortunately, repeating the 13 CO observations toward 
other pointings with high ratios would have required very long integrations. 
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(a) (b) 




Fig. 6. — A small, constant density, constant excitation spherical clump model, showing the 
geometry of the r = 1 surfaces for the CO (2 — > 1) and (1 — > 0) transitions for an observer 
located (a) to the far right of the clump let, and (b) in front of the clumplet. The very 
lightly shaded area surrounded by the dashed line is optically thin in both transitions and 
contributes a small (albeit not negligible) amount to the brightness of the clump in either 
transition. The darker shaded areas indicate the regions where first the (2 — >• 1) and then 
the (1 — > 0) transitions accrue enough column density to become optically thick. Most CO 
photons arise from these surfaces. Small, warm clumps fill more of the beam in the higher, 
more opaque CO transitions. 
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Fig. 7. — CO (2 — > 1)/(1 — > 0) ratio for clumps of different central (2 — > f) opacity r^, 
similar those of Figure 6, plotted against gas kinetic temperature. Note tha the overall 
opacity of an spherial, constant density clump is 33% higher than t c , which is the clump 
radius in opacity units. The calculations use LTE, assuming that the volume density is high 
enough to thermalize the lower CO transitions (n > 10 4 cm -3 ), and include the effect of the 
optically thin portion of the clump. The shaded areas represent the observational constraints 
for the region south of N83 where high ratios were found, as discussed in the main text ( 12 CO 
(2 -> 1)/(1 -»• 0)w 2, r 2i ~ f .5 - 4.4). 



